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ABSTRACT 

We show that the expansion of ejecta in Tycho's supernova remnant (SNR) is consis- 
tent with a spherically symmetric shell, based on Suzaku measurements of the Doppler 
broadened X-ray emission lines. All the strong Ka line emission show broader widths at 
the center than at the rim, while the centroid energies are constant across the remnant 
(except for Ca). This is the pattern expected for Doppler broadening due to expansion 
of the SNR ejecta in a spherical shell. To determine the expansion velocities of the 
ejecta, we applied a model for each emission line feature having two Gaussian compo- 
nents separately representing red- and blue-shifted gas, and inferred the Doppler velocity 
difference between these two components directly from the fitted centroid energy dif- 
ference. Taking into account the effect of projecting a three-dimensional shell to the 
plane of the detector, we derived average spherical expansion velocities independently 
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for the Ko emission of Si, S, Ar, and Fe, and K/3 of Si. We found that the expansion 
velocities of Si, S, and Ar ejecta of 4700 ± 100 km s _1 are distinctly higher than that 
obtained from Fe Ka emission, 4000 ±300 km s^ 1 , which is consistent with segregation 
of the Fe in the inner ejecta. Combining the observed ejecta velocities with the ejecta 
proper-motion measurements by Chandra, we derived a distance to the Tycho's SNR of 
4 ± 1 kpc. 

Subject headings: supernova remnants — supernovae: general — X-rays: individual 
(Tycho's SNR) 



1. Introduction 



Type la supernovae (SNe), the thermonuclear explosions of accreting C+O white dwarfs, are 
important for their role as standard candles to measure the distance to their host galaxies, and also 
as the main sources for the production of the Fe-group elements. However, our understanding of the 
physical mechanism of the explosion, as well as the details of the progenitor systems, is incomplete. 
The kinetic energy distribution of the ejecta holds important information about the propagation 
of the supernova burning front and the resulting nucleosynthesis. In sufficiently young supernova 
remnants (SNRs), this information is retained to some degree and can be probed by studying the 
X-ray emission of the shocked ejecta. 

Tycho's SNR is the remnant of the Galactic SN recorded by T ycho Brahe in 1572, and is 
secur ely classified as a Type la SN R based on the obser ved light curve (jBaaddll945l : iRuiz-Lapuente 
2004 ) and the light-echo spectrum (|Krause et al.ll2008l ). As the prototypical la SNR, Tycho's SNR 
has been targeted for study of the explosion mechanism. iBadenes et al.l (|2006l ) made a comparison 
between spatially integrated X-ray spectra of Tycho's SNR taken by XMM-Newton and Chandra 
and synthetic X-ray spectra based on several different Type la explosion models. They concluded 
that the observed spectra are well reproduced by a one- dimensional delayed detonation model 
with some degree of chemical stratification, with Fe-peak elements interior to intermediate-mass 
elements (IMEs, i.e., Si, S, Ar, Ca). 

The narrow band images of Tycho's SNR obtained by ASCA show that t he Fe-Ka emission 
is ac tually located i nterior to the emis sion lines of IMEs and Fe-L emission (jHwang &; Gotthelf 
19971 ). Furthermore. iHwang et al.l (|1998l ) found that the Fe-K emitting ejecta have a higher electron 
temperature and a lower ionization age than the IME ejecta. These results imply that the Fe-K 
emitting ejecta are situated more inside of the remnant and has b een heated by the r everse shock 
more rec ently than the other ejec ta. Observations with Chandra (IWarren et al.ll2005l ) and XMM- 
Newton ( Decourchelle et alJhoOll ) confirm the spectral segregation be tween the IME and Fe. The 
same situation has been reported fo r other la SNRs such as S N 1006 (jYamaguchi et al.1 12008 ) and 
LMC SNRs 0509-67.5 and 0519-69 (|Kosenko et alJl2008l . 1201(3 ) . 
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Thanks to the the spectral capability of Suzaku, iFuruzawa et al.l (|2009i ) discovered a significant 
broadening of Fe Kq line at the center of the Tycho remnant. They concluded that the shell of 
Fe Ko emitting ejecta is expanding with a line-of-sight velocity of 2800—3350 km s _1 . They 
also inferred line broadening for He-like Si and S Ko emission, but did not compute quantitative 
velocities. In this paper, we investigate the three-dimensional structure of the ejecta in Tycho's 
SNR and determine the expansion velocities. 



2. Observations and Data Reduction 

Suzaku observations of Tycho's SNR and an off-source background were carried out on 2006 
June 26-29 and 29-30, respectively, as a part of the Scientific Working Group observing time. The 
pointing position for the SNR was (RA, Dec) = (00 h 25 m 20 s , 64°08'18"); that for the background 
position was (00 h 36 m 54 s , 64°17'42"), 1?27 offset from the SNR along the Galactic plane. 

Suzaku carries two active instruments: four X-ray Imaging Spectrometers (XIS: Koyama et 
al. 2007) placed at the focal planes of four X-ray Telescopes (XRT: Serlemitsos et al. 2007) and 
a non-imaging Hard X-ray Detector (HXD: Takahashi et al. 2007). Each XIS has a 17/8x17/8 
field of view with a half-power diameter for the XRT of ~ 2'. One of the four XIS sensors (XIS 
1) is a back-illuminated CCD with high sensitivity at < 1 keV, while the others (XIS 0, 2, and 3) 
are front-illuminated (FI) CCDs with high efficiency and low background at > 5 keV. All of the 
XIS sensors were operated in the normal full-frame clocking mode without any spaced-raw charge 
injection (SCI: Uchiyama et al. 2009) for both the SNR and background observation^. 



The HXD data sanalysis were already reported by iTamagawa et al.1 (120091 ) , so we present here 
the detailed analysis of the XIS data. We focus on only the FI sensors (XIS 0, 2, 3), because these 
have better calibration of the energy gain near the Si-edge. We reprocessed the revision 2.0.6.13 
data products using the xispi software (HEAS0FT version 6.5) with the version 20080825 makepi 
file. After the reprocessing, the data were cleaned with standard screening criteria for the cleaned 
event data of the XIS. The effective exposures of the SNR and background were 101 ks and 51 ks, 
respectively. 

For the spectral fitting, we used response matrices created with xisrmfgen software version 
2007-05-14 using the version 20080311 rmf param file. We estimated the accuracy of the XIS energy 
scale for our own data and response functions using the line centroid energies of the 55 Fe calibration 
sources at two of the four corners of each XIS chip. We fitted Mn Ka lines from 55 Fe with Gaussian 
models and found that the line centroids of XIS 0, 2, and 3 were all 0.1—0.14% higher than 



1 Tycho's SNR was observed again by Suzaku on 2008 August 4-8 and 11-12 as one of the large proposal programs 
of the third announcement of opportunity observing cycle, this time with SCI on. However, we report only the 
analysis of the 2006 data in this paper becau se the current calibration of the response for the SCI-off data is better 
than for the SCI-on data (jOzawa et alJ 2009). 
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the expected energ y of 5.985 keV. Thi s is somewhat better than the systematic uncertainty of 
±0.2% reported by lOzawa et al.l (|2009i ). We also checked the width of the Mn Ka line is fully 
consistent with the expected instrumental broadening, even with all three FI sensors merged. The 
response files were generated by the xissimarf gen software, assuming a homogeneous disk-like 
radiation source with a radius of 4' corresponding to the azimuthally averaged radius of Tycho's 
SNR measured from the Chandra image. 



Overall Features 



A three-color XIS FI image of the SNR is shown in Figure [TJ Red, green, and blue colors 
correspond to the narrow energy bands of He-like Si Ka (1.7—2 keV), Fe Ka (6.2—6.7 keV), and 
the hard continuum band (7—13 keV). The images of the Si and Fe Ka emission, which represent 
the ejecta distribution, are brighter in north and fainter in south, while the hard band image 
is brightest in the sou thwest. These agree with the trends seen in more detaile d Chandra and 



XMM-Newton images (jHwang et al 



2002 



Warren et al. 



2005 



Badenes et alJl200fl ). 



The background-subtracted XIS spectrum of the entire SNR is shown in Figure [2j The source 
spectrum was taken from a circle with a radius of 4^65, with data from all three FI sensors merged 
to improve photon statistics. The background spectrum was extracted from a 7'-radius circular 
region on the off-source observation. In the spectrum, we can identify prominent K-shell emission 
line features of the He- and also some H-like ions of Mg, Si, S, Ar, and Ca, as well as emission from 
less ionized Cr, Mn and Fe (jTamagawa et al.ll2009i ). 



4. Spatially-Resolved Spectra 



Figure 4 in iFuruzawa et al.1 (120091 ) clearly shows that the Ka line blends of Si, S, and Fe are 
broadened at the center of the remnant, relative to the rim. However, their quantitative study 
was limited to the Fe Ka emission. Therefore, we extend their analysis to the lines of IMEs. We 
divided Tycho's SNR into 4 radial regions as shown in Figure [TJ with a circle at the center and 
three-quarter circular rings numbered from 1 (inner) to 4 (outer). The radius of the inner region 1 
was l'Al, and the thickness of the surrounding outer regions was l'.Q8. We excluded the southeast 
quadrant (60° — 150° where the angles increase counterclockwise from north ) because of the presence 
of irregular ejecta clumps (IVancura et al.lll995l : iDecourchelle et al.ll200ll ). We adopted the same 



background spectrum as used by Since the spectrum at energies below 1.7 keV is dominated by 
the Fe L- shell emission, where the atomic physics is complicated, we focused only on the energy 
band above 1.7 keV in the following sections. For simplicity, we divided the spectrum into two 
energy bands: the 1.7—5 keV band for the lines of IMEs ( §4.ip . and the 5—8 keV band for the Fe-K 
lines ( jgSp . 
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4.1. 1.7-5 keV Spectra 

To examine the radial changes of emission features, we fitted spectra taken from each re- 
gion with a phenomenological model featuring an absorbed power-law for the continuum emission 
plus twenty Gaussian components for the line emission. The column density for the fore ground 



interstellar gas was set to be Nn = 7 x 10 21 cm~ 2 following ICassam-Chenai et al.l ( 20071 ). The 
Gaussians represented transitions for the elements Si, S, Ar, and Ca for He-like and H-like ions: 
Hea (n = 2 — > n = 1), He/3 (ls3p — > Is 2 ), He7 (IsAp — > Is 2 ) in the He-like ions, and Lya (2p — > Is), 
Ly/3 (3p — > Is) in the H-like ions. We modeled the Hea transitions by a single Gaussian, which 
actually contains a blend of the forbidden, inter-combination, and resonance lines. For example, 
in Si Hea blend, the rest energies of the constituent lines are 1839.4 eV (forbidden), 1853.7 eV 
(inter-combinatio n), and 1864.9 eV (reson ance), with the blending giving an extra effective width 



of about 15 eV (jHwang fc Gotthelj 119971 ). We also represented each blend of He/3 and its sur- 
rounding satellite lines by a single Gaussian model. We note that the effective width for Si He/3 
is comparable to that of Hea, based on the non-equilibrium ionization (NEI) model version 1.1 
in XSPEC. The continuum spectra must be a combination of thermal and non-thermal emission, 
but these contributions are difficult to estimate separately (e.g., Tamagawa et al. 2009). We thus 
checked that the subsequent results for the Gaussian components do not change significantly even 
if we use a bremsstrahlung model for the continuum instead of a power-law model. 

The line centroid, width, and intensity of the prominent Hea blends of Si, S, Ar, and Ca, as 
well as the He/3 blends of Si and S, were fitted freely. The other emission lines are too weak to 
constrain the Gaussian parameters. Thus, we fixed the energy difference between the prominent 
line features above and weaker line features of the same element; for example, the energy difference 
between S Hea and S Lya was fixed to 178 eV. Here we assumed electron te mperature kT P of 



1 keV and ionization timescale n e t of ~ 10 11 cm~ 3 s, typical in Tycho's SNR ( Hwang et al.lll998l ). 
because the centroid energy of the Hea blend depends on kT e and n e t; and used the Astrophysical 
Plasma Emission Database (APED: Smith et al. 2001). The line widths were linked to each other 
for lines of the same element (e.g., S Hea, S Lya, S He/3, S He7, and S Ly/3). The intensities of the 
prominent He- like emission above and the relatively strong Lya lines of Si and S, were freely fitted. 
For lines other than those, we fixed the intensity ratios: we took He7/He/3=0.3 and Ly/3/Lya=0.1 
for the plasma in Tycho's SNR (kT e =l keV and n e t=10 11 cm" 3 s, as noted above). These intensity 
ratios do not vary much with n e t and vary only 10% over a decade in kT e ; hence the uncertainties 
here are negligibly small. We also fixed the ratios of Lya/Hea and He/3/Hea for Ar and Ca to be 
the same as these of S. 

With these models and assumptions, we fitted all the four spectra. The reduced x 2 values and 
degrees of freedom (dof) in regions 1, 2, 3, and 4 were 1.37 (791), 1.42 (724), 1.74 (858), and 1.56 
(846), respectively. The large reduced x 2 values are attributed to the high statistics of spectra and 
our approximations above for the fitting model. In addition, there might be some contaminating 
emission from the less ionized ions (e.g., Li-like). We thus conclude that the best-fit models do 
reproduce the spectra sufficiently well. Individual components of the best-fit models for regions 1 
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and 4 are shown in Figure[3] (a-1) and (b-1), respectively. The best-fit parameters of all the fits are 
listed in Table [TJ 



4.2. 5 8 keV Spectra 

We also fitted the 5 — 8 keV band spectra with a power-law for the continuum, plus three 
Gaussian components for Cr and Fe Ka, and Fe K/3 line blends. Since Cr Ka and Fe K/3 emission 
are very weak, these widths were tied to that of Fe Ka. The other parameters were allowed to vary 
freely. The reduced x 2 s and dof for regions 1, 2, 3 and 4 were 1.25 (115), 1.22 (214), 1.07 (335), 
and 0.86 (298), respectively. All the fits for the four spectra from regions 1—4 were acceptable. 
The 5—8 keV spectra from regions 1 and 4 with the best-fit models are shown in Figure [3] (a-2) 
and (b-2), respectively. The best-fit parameters are listed in Table [TJ 



5. Radial Line Profiles 

Figure [4] shows the centroid energies, widths, and intensities of the Ka line blends of Si, S, 
Ar, Ca and Fe, obtained in the previous section. We found that the centroid energies of all the 
elements except for Ca are constant with respect to radius within the sys tematic un c ertain ties of 



±0.2%, indicated in each left box. This systematic error was estimated by lOta et al.l (120071 ) using 
the Suzaku observations of extended sources with Fe Ka emission (Cygnus Loop and Sagittarius C), 
and represents the gain variation on the same CCD chip at the 90% confidence level. In contrast to 
the constancy of the centroid energies, the widths of all the emission features decrease significantly 
from the center to the rim. The difference between the Si and Fe widths of regions 1 and 4 is 
10.5 ± 0.3 eV and 33 ± 6 eV, respectively. 

The radial intensity profile of IMEs emission are similar to each other with a peak radius of 
3' — 4', while Fe-Ka emission has a radius of ~ 3' which is somewhat smaller than that of the IMEs. 
These trends agree with the previous observations (e.g., Hwang & Gotthelf 1997). 

There is some indication of spectral variations within our source regions in Figure [lj so we 
investigated the azimuthal variations of the line centroids and widths by dividing regions 2, 3, 
and 4 evenly into three azimuthal sectors. We applied the same fitting procedure to the spectrum 
accumulated from each divided region. There is some variation amongst the azimuthally divided 
sectors, with the azimuthal variations of line centroids being ±3 eV for Si, ±5 eV for S, ±8 eV for 
Ar, ±20 eV for Ca, and ±10 eV for Fe. These values are much smaller than the corresponding 
line widths; however, in all the radial sectors, the widths show the same trends of being wider in 
the inside compared to the rim, and the centroids being more constant. The qualitative trends 
appear to be the same as for the azimuthally averaged spectral regions, therefore we proceed with 
our azimuthally averaged analysis. 
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In the following subsections, we consider two simple interpretat ions of the radial properties of 
centroid energies and widths as discussed by iFuruzawa et al.l (|2009l ): (a) plasma characterized by 
multiple ionization ages, and (b) Doppler shifts resulting from expansion in a shell. 



5.1. Plasma with Multiple Ionization Ages 

It is expected that the plasma in the SNR will have multiple ionization ages. Since the 
reverse shock propagates from the outside of the remnant inward, the shocked ejecta will have 
different n e t between the outside (near the contact discontinuity) and the inside (where the reverse 
shock has just passed through), with the more recently shocked ejecta generally having a lower 
n P t and a different kT e from t he outer ejecta, with details depending on the ejecta density profile 



(jDwarkadas fc Chevalierlll998l ). 



Plasma with a range of n e t values easily makes the observed line blend broaden at the center 
of the remnant, since line centroids strongly depend on n e t. For example, based on the NEI model 
of version 1.1, the centroid of the Ka line blend of S at kT e =l keV varies from 2.41 keV to 2.45 
keV depending on the n e t values of 10 10 — 10 11 cm -3 s, which may result in a line broadening 
of ~40eV. However, our results show that the centroid energies are constant across the remnant 
within 0.2%. This implies that the plasma emitting each line or blend has the same average n e t all 
over the SNR. Therefore, a multiple n e t plasma does not provide a fully satisfying explanation of 
our observation. 



5.2. Doppler Shift by the Shell Expansion 

Another case we consider is an expanding shell of ejecta. The emission lines from retreating and 
approaching gas are red- and blue- shifted, respectively. Thus, the line emission at the projected 
center of the remnant is expected to be broadened. At the rim, on the other hand, a narrow line 
would be observed, since there is only a small line-of-sight component of the velocity. In addition, 
if the expansion is spherically symmetric (and the plasma is reasonably uniform), the centroid 
energies will be constant with radius. We found that the observed line profiles are well reproduced 
by the Doppler broadening of the spherically symmetric shell expansion. 



5.3. Properties of Ca He-like Ka line 

As shown in Figure HI the Ca Hea blend is the only feature where the centroid energy is 
seen to gradually increases from the center to the rim. Possibly, this suggests that the Ca ejecta 
might be in multiple ionization ages unlike the other ejecta. However, the Ca Hea blend is heavily 
contaminated by Ar He/3, and thus the inferred properties of the Ca emission depends rather 
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strongly on the assumed intensity ratio, He/3/Hea of Ar. A future observation with higher energy 
resolution is required to reliably investigate the accurate properties of Ca Hea blend. 



6. Velocities of Ejecta 



In this section, we quantitatively derive the expansion velocity with the Ka emission of Si, S, 
Ar, and Fe, and K/3 of Si, assuming that the ejecta are expanding in a spherical shell. 



6.1. Method of Velocity Determination 

In a spherically symmetric shell expansion, a broadened line should have both red- and blue- 
shifted components. If we apply such a two-Gaussian model (with red- and blue-shifted compo- 
nents) to a single broadened line, the line energy shift (i.e., relative to the rest frame) gives the 
Doppler velocity as 

\E bs,i - Ep\ _ V±j 

hiQ c 

where -E bs,i and v±i are the observed centroid energy of the red- or blue-shifted line and the 
line-of-sight velocity in each region i (i = 1 — 4), respectively, and Eq is the line energy at the 
rest frame. We introduce here the energy shift in each region i, SE^ = l-Epbs,* — Eq\. Assuming 
front-back velocity symmetry, the parameter we measure is 2 x 5Ei, corresponding to the energy 
difference between the red- and blue-shifted lines. 

Once we obtain v±i, we can convert it to the expansion velocity v exp by considering two 
important effects: the projection of the three-dimensional shell onto the plane of the detector, and 
the limited spatial resolution of the Suzaku XRT. In an appendix, we calculate the factors which 
represent what percentage of the u ex p would be observed as v±i in each region i. In other words, 
v exp can be expressed in terms of and the observable parameter v±i as 

Vexp = £7 ■ ( 2 ) 

The calculated Cj factors are summarized in Table [2j 



6.2. Spectral Fitting 

As a first step, we derived the expansion velocities by fitting the spectra of regions 1 and 4 at 
the same time, because the difference between v±i and u_i_4 is expected to be the largest. 

To test for spherically symmetric expansion, we therefore fitted the 1.7—5 keV spectra with the 
model of an absorbed power-law for the continuum and Gaussian components for emission lines. 
The column density for the foreground interstellar gas was set as in §4.1 1 We applied Gaussian 



- 9 - 



models including red- and blue-shifted components, for each line feature. We left the centroid 
energies of the both lines free, while the widths and intensities were set to be equal to each other. 
Because we fitted the spectra of regions 1 and 4 at the same time, we actually applied two pairs 
of the red- and blue-shifted lines (four Gaussians in total) for each line feature. According to 
Equation [H the ratio of dEi/SE^ equals the ratio of w_li/Vi_4- Since v exp should be the same for 
regions 1 and 4, Equation [2] gives 5Ei/5E^ = C1/C4. In our fitting, we therefore fixed the ratio of 
5E 1 /SE 4: = C1/C4. As a result, 

^expj the width, and intensity were the only three parameters that 
were allowed to vary freely for each set of four Gaussians. 

We employed twenty sets of four Gaussians for the observed line features: Hea, Lya, He/3, 
He7, Ly/3 of Si, S, Ar, and Ca. We derived the v exp independently for the prominent Hea blends 
of Si, S, and Ar, and He/3 of Si. The energy separation between these and other line features were 
fixed as the model in §4.1} for example, the energy difference between Ar Hea and Ar He/3 was 
fixed to 209 eV (see Table [JJ. The sole exception is S Lya, whose energy was free in §4.1 1 here, we 
fixed the energy shift from S He/3 to 181 eV to compensate for limited photon statistics. A single 
line width was fitted for lines of the same element (e.g., for S: Hea, Lya, He/3, He7, and Ly/3). The 
exception was Si Hea, whose width was fitted independently from that of other Si lines, due to 
uncertainties in the response near the Si-edge. As noted in £14.11 the width of the Gaussian applied 
to the Hea blend includes a contribution from the blending of the constituent triplet lines (~ 15 
eV). This broadening varies within a few eV depending on the plasma condition, but we checked 
that the measurement of the line centroid is little affected by the uncertainty of the width. The 
intensity ratios of weak lines to prominent lines (e.g., S Lya/Hea) were fixed by assuming a plasma 
with kT e r-j 1 keV and n e i~ 10 11 cm ~ 3 s, as in §4.11 The exception here is that S Lya/Hea was 
fixed to the value obtained from the fits in §4.11 (Table [T]) . 

We also fitted the 5 — 8 keV spectra of regions 1 and 4 using the same procedure as above. 
The model consisted of a power-law continuum and red- and blue-shifted Fe Ka lines. We also 
included the Cr Ka and Fe K/3 blends, but as single Gaussians due to limited photon statistics. We 
fixed the ratio of bE\jbE^ to C1/C4 following Table El The widths and intensities of the red- and 
blue-shifted lines in regions 1 and 4 were all linked to each other. Figure [5] shows the individual 
components of the best-fit models. The reduced x 2 s and dof were 1.47 (1430) for 1.7—5 keV and 
0.96 (415) for 5—8 keV band spectra. 

As a second step, we have also fitted all the spectra of regions 1, 2, 3, and 4 at the same time 
in exactly the same manner just described. Here our aim was to verify that the model was valid for 
all four regions, and also to improve the photon statistics. The ratios of energy difference among 
regions 1, 2, 3 and 4 were fixed as 5E\ : 5E2 : 5E% : 5E4 = C\ : C2 ■ C3 : C4. The parameters 
obtained from the fits are consistent with those obtained for regions 1 and 4, and within the 
systematical errors. The reduced x 2 s an d dof were 1.55 (3028) and 1.06 (970) for 1.7—5 keV and 
5—8 keV band spectra, respectively. The best fit parameters are summarized in Table El 
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6.3. Expansion Velocities 



We translate 5Ei into v±i using Equation [fl and then convert to v exp using Equation [2] and C, 
(Table [2]). The derived velocities are summarized in Table We note that the ? w p of the Fe Kg 



Furuzawa et al 



emiss ion in Table [3] does not correspond to the line-of-sight velocity measured in 
(|2009l ). because the v exp represents the expansion velocity in which the projection effect and PSF 
effect have properly been taken into a ccount. Meanwh i le, th e v±\ is the value should correspond to 
the line-of-sight velocity measured in iFuruzawa et al. ( 200S), 30 4 km s — 1 . Our measurement 



of v±i and the line-of-sight velocity measured in lFuruzawa et al.l (|2009l ) agree well each other. The 



J exp 



derived from Si Hea and He/3 are consistent with each other. Therefore, we conclude that the 
systematic uncertainty of the Si-edge in the response is small enough to determine the velocity. 
The expansion velocity obtained from the Hea blends of Si, S, and Ar are clearly higher than that 
obtained from the Fe Ka blends. 



7. Discussion 

Based on the analysis of spatially resolved spectra, we have demonstrated that (1) the ejecta 
shell is expanding in a generally spherical and symmetric manner, and (2) the expansion velocities 
of the IME ejecta are significantly higher than that of the Fe Ka emitting ejecta. These results 
then allow us to discuss the distance to Tycho's SNR, the ejecta segregation, and the reverse shock 
velocity. 



7.1. Distance to Tycho's SNR 



Using Ch a ndra high-resolution images of Tycho's SNR obtained in 2000, 2003, and 2007, 
Katsuda et al.l (|2009l ) measured the expansion rates of both the forward shock and the reverse 
shocked ejecta. They found tha t the mean pro per-motion of the Si-rich layer is /i ~ / .'25 yr _1 . 
This is consistent with the result iHughesI (120001 ) had derived using ROSA T data. Combining with 
our expansion velocities of 4700 ± 100 km s" 1 (Table EJ), we obtained a range of the distance to 
the SNR of D = (4.0 ± 1.0)(u/4700 km s -1 ) (fi/tf! 25 yr -1 ) -1 kpc. The relatively large uncertainty 
is mainly due to the azimuthal variation of the proper motion of the ejecta. Our result is the first 
estimate of the distance to Tycho derived solely from X-ray observations. 

The distances inferred by previous studies is shown in Figure [6j The distance of around 2—3 
kpc estimated by modeling the observed Ha line spectra (gr een in Figure [6|) has been most widely 
cited thus far, but is model-dependent. iKrause et al.1 (120081 ) recently derived a larger distance of 
3.8jli'i kpc, based on the SN peak luminosity est imated by the observ ed optical light-echo spectrum. 
Our estimate is consistent with the result from lKrause et aO ( 2008 ). 
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7.2. Ejecta Segregation 



The Fe-K radial profile peaks at 
and Fe-L emission of 



180" which is a distinctly smalle r radius than the IME 



200" (Decourchelle et al 



200f 



Warren et al.ll2005l ). In addition, the Fe-K 
emitting ejecta has about 100 times lower n e t than the IMEs and Fe-L emitter (10 11 cm s _1 : Hwang 
et al. 1998). These results suggest that the plasma emitting Fe-K has been heated by the reverse 
shock more recently than the other materials. Since the reverse shock propagates from the outside 
toward the center, the Fe-K and Si emission are tracers of the interior and exterior material in the 
SNR, respectively. 

Our results show that the Si ejecta also have higher f exp (4700 ± 100 km s _1 ) compared to the 
Fe Kq emitting ejecta (4000 ±300 km s _1 ). These velocity measurements add to the morphological 
and spectral evidence that the Fe-K emitting ejecta are segregated interior to the Si ejecta. The 
expansion rate ratio of Fe Ko emission to Si ejecta is derived to be (4000/ 180)/ (4700/200) ~ 0.93. 



7.3. Reverse Shock Velocity 



The angular radius of the reverse shock of 183" estimated by I Warren et al.l (|2005l ) corresponds 
to r rs — 3.5 i 0.9 pc at a distance of 4.0 ± 1.0 kpc. Given the age of 434 yr, the velocity of 
unshocked ejecta f U n,ej is then estimated to be 8000 ± 2000 km s _1 at a radius of 3.5 pc. The 
post-shock velocity f s h,ej can be approximated as the obtained expansion velocity of Fe, 4000 ± 300 
km s . Under the assumption of strong shock, conservation laws for mass, momentum, and energy 
across the shock front yield the following Rankine-Hugoniot relation, 



un i e J 



7 + 1 



7 



1 



(Vsh,, 



(3) 



where t> rs and 7 ar e the reverse s hock velocity in the observer's frame and specific heat ratio, 
respectively. Since I Warren et al.l (120051 ) argued that the reverse shock of Tycho's SNR is not 
strongly accelerating cosmic-rays (unlike the forward shock), we here take 7 = 5/3 as appropriate 
for a non-relativistic ideal gas. Equation [3] then gives v IS = 2700 ± 800 km s , and hence the 



upstream velocity in the shock-rest frame is obtained to be v T 



un i e j 



5300 ± 2100 km s" 



This value allows us to estimate the explosion energy of Tycho's SN by comparing with a 
self-similar model for the evolution of you ng SNRs interacting with a uniform density ISM. In 
Figure 2d of bwarkadas fc Chevalien (Il998h . v rs is expected to be 3.4 x 10 3 E^ 2 (M e] /1AM Q )-^ 2 
km s _1 throughout most of the SNR's evolution, for either r~ 7 power-law and exponential density 
distributions of the ejecta, where E§\ is the explosion energy in units of 10 51 ergs. Thus, assuming 
an ejecta mass of 1.4Mq, the explosion energy of Tycho's SN is derived to be (2.5 =L 2.0) x 10 51 
ergs. Given that the error is large, the derived energy range includes the standard value for normal 
Type la SNe. 



From the Rankine-Hugoniot relations, the post-shock temperature for each particle species a 
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is given as 



kT n 



2(7 ~ 1) -2 _3_ _ 2 

try + l\2 maV ™ - ^""WsJ 



(4) 



where m a is the particle mass. Given v rs = 5300 ± 2100 km s _1 , the temperature of the shocked 
Fe is expected to be kTp c = 3 ± 2 MeV. Such a high temperature would result in a line broadening 
of AE = 70 ± 30 eV due to the thermal Doppler effect. We measured the width of the Fe-K blend 
of 55±3 eV (Table [3]), where an intrinsic line broadening of 20 — 30 eV is expected for the plasma 
in th e range of kT e > 1.6 keV and 10 < log n e t cm~ 3 s < 10.7 (jHwang et al.lll998l : iFuruzawa et al 



2009). The additional width of ~50 eV which cannot be explained by the intrinsic line broadening 



might be explained by the thermal Doppler broadening of Fe ions. Future missions with high energy 
resolution, like Astro-H, will help us to accurately measure ion temperatures and study the heating 
process at the reverse shock of young SNRs in detail. 



8. Summary 



We analyzed the X-ray spectra of Tycho's SNR using the Suzaku 100 ks observation. We 
summarize our results below: 



1. We obtained the radial dependence of centroid energies, widths, and intensities for the Ka 
emission of Si, S, Ar, Ca, and Fe for annular regions covering nearly three- fourths of the 
azimuth around the SNR. The centroid energies of all line blends except Ca are constant 
with radius, while the widths significantly become narrower from the center to the rim. We 
found that the observed line properties are well explained by a spherically symmetric shell 
expansion of the ejecta. 

2. We derived the expansion velocities of ejecta with Doppler broadened Ka emission of Si, S, Ar 
and Fe, and K/3 of Si independently. We fitted each broadened line feature with two Gaussians 
representing red- and blue-shifted gas, and obtained the expansion velocity from the centroid 
energy separation between two lines. We found that expansion velocities measured for Ka 
emission of Si, 4700 ±100 km s -1 , are clearly higher than that measured for Fe Ka, 4000 ±300 
km s^ 1 . 

3. Combining the obtained expansion velocity with the proper-motion measurement by Chandra, 
the distance to Tycho's SNR is estimated to be 4 ± 1 kpc. 

To determine the true expansion velocity f exp , we have to calculate the correction factors Cj 
introduced in Equation [2j C{ describes the projection effect of the three-dimensional shell onto the 
detector plane, including the photon smearing effects caused by the limited spatial resolution of 
the Suzaku XRT. In this appendix, we calculate the correction factors Cj. 



-13- 



First we define the regions that we use in the following calculations. As shown in Figure [7J 
we divide the entire SNR into ten "SKY" regions. The region enclosed by the innermost circle is 
named SKY 1, the next two annular regions are SKY 2 and 3, and the outer six annular regions, 
in which a quarter region is excluded, are named SKY 4 (inner) to SKY 9 (outer). The radii of the 
SKY regions are 30" for SKY 1, and are incremented by 30" toward the outer edge of the SNR. 
The quarter region in the southeast is named SKY 0. The regions 1 — 4 are the detector regions 
used in Sj4] (Figured]). The regions 1, 2, 3, and 4 roughly correspond to the SKY 1 — 3, 4 — 5, 6 — 7, 
and 8 — 9, respectively. 



A. Projection Effects 

Consider a thin and spherically symmetric shell with radius r s ^. The shell is expanding with 
velocity v exp . The polar angle 9 is defined as the angle measured from the horizontal line which 
intersects the center of the spherical shell, and the azimuth angle <p represents an angle around the 
horizontal line. The sky region j (j = 0, 1, . . ., 9) is described by using the range of 9 (0{ and 9 2 ) 
and 4>. The 9\ and 9 2 angles are given as 

9\ = arcsin (r 3 /r sh ) , 9% = arcsin (r 9 /r sh ) , (Al) 
9{ = arcsin (rj_i/r sh ) , 9 3 2 = arcsin (rj/r sh ) (j = 1,2, .. .,9), (A2) 

(A3) 

where the radius r& is given as 

r fc = fc-30" (Jfe = 0,l,...,9). (A4) 
We here consider only the half of the spherical shell (0 < 9 < 90 deg) in to evaluate 6E{. 

The projected velocities Vj of the SKY region j are found as the mean line-of-sight velocity 
v± = v exp cos 9. By integrating v± over the solid angle, the mean velocity Vj is described as 

v j = V ex p ' Aj , (A5) 

where 

cos 9 sin 9 d9 

1 WJ hi • (A6) 



sin 9 d9 

Since 9\ and 9\ are functions of r s h, Aj also depends on r^. We calculated Aj for the radius r s h at 
2 arcsec intervals. 



The actual ejecta shell has a finite width in radius r s h- iDecourchelle et al.l (|200ll ) and lWarren et al 



(|2005T ) reported r sh = 190" - 220" for Si, and r sh = 180" - 200" for Fe. We therefore find an average 
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of Aj over the given range of r s h as 



min 
sh 

—max 
'sh 



Aj = • (A7) 

'sh 



In Table HI we summarize the calculated A,-. 



B. Limited spatial resolution of XRT 

We estimate the smearing effect caused by the limited spatial resolution of the XRT. We use 
the xissim software, which simulates interactions of X-ray photons incid ent upon the XRT/XIS 



system using ray-tracing and Monte Carlo techniques (llshisaki et al.l 120071 ) 



First, we prepare the input image files for the simulation. Because the distributions of K shell 
emission for IMEs and Fe are slightly different (e.g., Hwang & Gotthelf 1997), we prepare two 
narrow-band Chandra images: the Si and S (1.8 — 1.9 plus 2.4 — 2.5 keV) band, and the Fe-K (6 — 7 
keV) band. Next, we divide each Chandra image into the SKY regions introduced in Figure and 
we take this divided image as an input to the xissim simulator. In the simulation, we fix the input 
photon energy to be monochromatic: 1.8 keV and 6.5 keV for Si+S and Fe, respectively. Finally, 
in the simulated images, we accumulate the number of photons in each detector region 1 — 4. We 
repeat this procedure for both energy bands (Si+S and Fe) and for all ten different SKY regions 
(SKY 0-9). 

Results of the simulations are summarized in Table [5l We introduce the smearing factor F{j , 
which represents the fraction of the number of photons that originate in a certain SKY region j 
( j = — 9) relative to to the number of photons detected in a certain detector region i (i = 1 — 4). 
For example, in the Fe K band, = 0.194 means that 19.4% of photons detected in region 1 
originate from SKY 3. 



C. Correction Factors 

We derive the correction factor Cj by multiplying the projection factor Aj (Table by the 
smearing factor Fij (Table [5]) as 

9 

Ci = J2 F ij A 3 (» = 1,2,3,4). (CI) 

3=0 

The derived values are summarized in Table [2j 
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Fig. 1. — Three color XIS FI (XIS 0+2+3) image of Tycho's SNR. Red, green, and blue correspond 
to the energy bands of Si Hea (1.7—2 keV), the Fe Ka (6.2—6.7 keV), and the hard band energy 
band (7—13 keV), respectively. The white circle and quarter sectors are the regions where we 
extract spectra. 
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Fig. 2. — Background-subtracted XIS FI (XIS 0+2+3) spectrum of the entire SNR. The important 
lines are marked. Hea, He/3, and Lya indicate emission lines of n = 2 — > n = 1 transition in He-like 
ion, ls3p — > Is 2 transition in He-like ion, and n = 2 — > n = 1 transition in H-like ion. 
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Fig. 3. — Background-subtracted 1.7—5 keV (left) and 5—8 keV (right) band spectra extracted 
from region 1 (a) and 4 (b) indicated in Figure [H The best-fit models are shown in solid lines and 
residuals are given in the bottom panels. 
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Fig. 4. — Radial dependence of the best-fit Gaussian parameters of Hea lines of Si, S, Ar, and Ca, 
and Fe Ka blend. The centroid energies (left), the 1 a widths (middle), and the intensities (right). 
The systematic errors from the intrachip gain variation at 90% confidence level are represented at 
the right corner in each left box. The statistical errors also represent 90% confidence level. 
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Fig. 6. — Estimates of the distance to Tycho's SNR. 




Fig. 7.- 



SKY regions defined in the plane of the sky. The underlying image was taken by Chandra. 
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Table 1. Best-fit Gaussian parameters 



Region 1 Region 2 

Line Ccntroid Width (la) Flux 5 Ccntroid Width (Ict) Flux a 

(kcV) (eV) (kcV) (cV) 

Si Hea 1.8661+j;;™'^ 31.0±0.2 115.7±0.3 ^866^^3^ 27 -°-o.2 140.3+JJ'^ 

Si Lya = He/9 - 0.178 = He/3 3.5±0.1 = He/3 - 0.178 = He/3 ^S^og 1 

Si He/3 2.1972^° °°^ 36±1 10.4+0.1 2.1962±° ) ;°°°J 30. 5t°' 6 13.26+gg 11 

Si He7 =Hc/3 + 0.111 = He/3 (He7/Hc/3=0.3) = He/3 + 0.111 = He/3 (He7/He/3=0.3) 

Si Ly/3 =Hc/3 + 0.193 = He/3 (Ly/3/Lyg=0.1) = He/3 + 0.193 = He/3 (Ly/3/Lyg=0.1) 

S Hea 2.4553±"£^£ 403^ 27.4+0.2 2.4549+0.0003 36.1+'^^ 36.3+0.2 

S Lya = He/3-0.263 = He/3 0.18+0.06 = He/3 - 0.263 = He/3 0.395j°, 4 ° 6 56 

S He/3 2.899+0.003 59±| 2.05+0.07 2.887+0.002 44l^ 2.58+0.05 

S He7 =Hc/3 + 0.149 = He/3 (Hc7/Hc/3=0.3) = He/3 + 0.149 = He/3 (He7/He/3=0.3) 

S Ly/3 = He/3 + 0.222 = He/3 (Ly/3/Lyg=0.1) = He/3 + 0.222 = He/3 (Ly/3/Lyg=0.1) 

Ar Heg 3.130+0.002 49^ 2.351^'^ 3.129+0.001 47^ 3.44l"°;°| 

Ar Lya = Hea + 0.209 = Hea b = Hea + 0.209 = Hea b 

Ar He/3 = Hea + 0.561 = Hea c = Hea + 0.561 = Hea c 

Ar He7 = Hea + 0.751 = Hea (Hc7/He/3=0.3) = Hea + 0.751 = Hea (He7/He/3=0.3) 

Ar Ly/3 = Hea + 0.812 = Hea (Ly/3/Lya=0.1) = Hea + 0.812 = Hea (Ly/3/Lya=0.1) 

CaHca 3.850+0.006 82±g 0.78+0.03 3.8571"°,; ^| 74+4 1.13+0.03 

Ca Lya = Hea + 0.220 = Hea b = Hea + 0.220 = Hea b 

Ca He/3 = Hea + 0.654 = Hea c = Hea + 0.654 = Hea c 

CaHe7 = Hea + 0.938 = Hea (Hc7/Hc/3=0.3) = Hea + 0.938 = Hea (He7/He/3=0.3) 

CaLy/3 = Hea + 0.979 = Hea (Ly/3/Lya=0.1) = Hea + 0.979 = Hea (Ly/3/Lya=0.1) 

Cr Ka 5.46+0.04 = Fe Ka 0.05+0.02 5.51+QQ5 = Fe Ka 0.05+0.02 

Fe Ka 6.445+0.004 86+4 1.00+0.03 6.4441"° °°l 71+3 1.37+0.03 

Fe K/3 7.14+0.05 = Fe Ka 0.05+0.02 7.11+0.02 = Fe Ka 0.08l°,;°,i 



Region 3 Region 4 

Line Ccntroid Width (la) Flux a Ccntroid Width (la) Flux a 

(kcV) (eV) (kcV) (eV) 

Si Hea O6504Z— 2^^ 064§F™ m^TTT lgl 2± o.3 

Si Lya = Hc/3 -0.178 = He/3 7.28±^ ?| = He/3 - 0.178 = He/3 7.5+0.1 

SiHc/3 2.19451°;°°°? 23.61°;? 17.691°;^ 2.1935l°;°°°| 21. 1±°J 18.0+0.2 

Si He 7 =Hc/3 + 0.111 = He/3 (He7/He/3=0.3) = He/3 + 0.111 = He/3 (Hc7/He/3=0.3) 

Si Ly/3 =Hc/3 + 0.193 = Hc/3 (Ly/3/Lya=0.1) = He/3 + 0.193 = He/3 (Ly/3/Lya=0.1) 

S Hea 2.4539+0.0002 29.9+°'j 49.6+°'^ 2.4535%>°°f 02 28.7+°;^ 50.7+°;^ 

S Lya =Hc/3 - 0.263 = He/3 0.75±° °!? = He/3 - 0.263 = He/3 0.761"°'°° 

SHc/3 2.8881°;"? 34l 2 3.54t°£ 2.891°,"°? 351°;^ 3.74±°£ 

S He7 = He/3 + 0.149 = He/3 (He7/He/3=0.3) = He/3 + 0.149 = He/3 (He7/He^=0.3) 

S Ly/3 = Hc/3 + 0.222 = Hc/3 (Ly/3/Lya=0.1) = He/3 + 0.222 = He/3 (Ly/3/Lya=0.1) 

Ar Hea 3^1273^009 4 - 84 ±o'.o6 3.127+0.01 38+3 5.111°;°,? 

Ar Lya = Hea + 0.209 = Hea b = Hea + 0.209 = Hea b 

Ar Hc/3 = Hea + 0.561 = Hea c = Hea + 0.561 = Hea c 

Ar He7 = Hea + 0.751 = Hea (Hc7/Hc/3=0.3) = Hea + 0.751 = Hea (Hc7/He/3=0.3) 

Ar Ly/3 = Hea + 0.812 = Hea (Ly/3/Lya=0.1) = Hea + 0.812 = Hea (Ly/3/Lya=0.1) 

Ca Hea 3.86121°,;","^ 62~f| 1.571^;°^ 3.868+0.03 49^ 1.52t°;°| 

Ca Lya = Hea + 0.220 = Hea b = Hea + 0.220 = Hea b 

Ca He/3 = Hea + 0.654 = Hea c = Hea + 0.654 = Hea c 

CaHe 7 = Hea + 0.938 = Hea (He7/Hc/3=0.3) = Hea + 0.938 = Hea (Hc7/He/3=0.3) 
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Table 1 — Continued 

CaLy/3 = Hea + 0.979 = Hea (Ly/3/Lya=0.1) = Hea + 0.979 = Hea (Ly/3/Lya=0.1) 

Cr Ka 5.48±0.03 = Fe Ka 0.06±0.02 5 - 44 to'.o4 = Fe Ka 0.05±0.03 

Fe Ka 6.448±+0.002 58+3 1.47±0.03 6 - 450 ±o'.oo4 53 -t L09 iao3 

Fe K/3 7.12±0.03 = Fe Ka 0.05^ 02 7.13+0.09 = Fe Ka 0.03+0.02 



Note. — Errors indicate the 90% confidence limits. 
a Units in xlO -5 photons cm -2 s — 1 arcmin - 2 . 
b ya/Hea is linked to that of S. 
c He,9/He« is linked to that of S. 



Table 2: Correction factors Cj denned in Equation [2] 









IMEs 


Fe 


Rej 


don 


1 


0.748 


0.712 


Rej 


;ion 


2 


0.578 


0.530 


Re^ 


lion 


3 


0.315 


0.262 


Rej 


don 


4 


0.170 


0.120 



Table 3: Best-fit parameters of the double Gaussian model 





Width 


-E-red 




25E Y 


v±i 


^cxp 




(eV) 


(keV) 


(keV) 


(eV) 


(km s" 1 ) 


(km s _1 ) 


Si Hea 
Si He/3 


20-16^ 


1 SSS'3+0.0002 

9 lfiSK+0.0012 
Z.1DOO_ 0007 


1.8443 ± 0.0002 

9 91 Q/l +0.0006 
Z - Zi3 ^-0.0014 


44.0t°| 
51^ 


3540±20 
3480i?° 


4730±|g 
4700±100 


S Hea 


26 - 6+ 03 


9 /197n+0-0005 
z '^ z '"-0.0004 


9 /IS/I1+ 0003 
2.4841_ 00005 


57.1±0.6 


3490±40 


4660±50 


Ar Hea 


34. 6 ; ; 4 7 


3.089±0.002 


3.164 ±+0.002 


74±3 


3600±100 


4800±200 


Fe Ka 


55±3 


a r; n7 +0.005 
D - ou ' -0.007 


fi QOQ+0.007 
O.oOo_ 005 


124±8 


2900±200 


4000±300 



Note. — Errors indicate the 90% confidence limits. 
"Best fit values of region 1. 

Table 4: Calculated Aj described in Equation IA7I 





SKY0 


SKYf 


SKY2 


SKY3 


SKY4 


SKY5 


SKY6 


SKY7 


SKY8 


SKY9 


Si Shell* 


0.455 


0.995 


0.973 


0.927 


0.854 


0.744 


0.573 


0.138 








Fe Shell 1 " 


0.447 


0.994 


0.968 


0.914 


0.827 


0.693 


0.450 












* r sh = 190" - 220' 
f r sh = 180" - 200' 
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Table 5: Results of the simulation: fraction of the photons detected in REG i to those originated 
from SKY j 

Si and S 

SKY SKY 1 SKY 2 SKY 3 SKY 4 SKY 5 SKY 6 SKY 7 SKY 8 SKY 9 

0.100 0.068 0.190 0.207 0.152 0.098 0.081 0.071 0.030 0.003 

0.047 0.010 0.034 0.087 0.171 0.203 0.195 0.175 0.071 0.006 

0.033 0.002 0.007 0.016 0.038 0.089 0.220 0.383 0.196 0.017 

0.031 0.001 0.003 0.007 0.017 0.037 0.104 0.325 0.393 0.082 

Fc 

SKY SKY 1 SKY 2 SKY 3 SKY 4 SKY 5 SKY 6 SKY 7 SKY 8 SKY 9 ~ 

0.105 0.064 0.178 0.194 0.154 0.113 0.099 0.062 0.024 0.006 

0.050 0.009 0.032 0.079 0.164 0.223 0.235 0.143 0.054 0.013 

0.037 0.002 0.006 0.016 0.039 0.100 0.271 0.335 0.160 0.036 

0.036 0.001 0.003 0.006 0.016 0.040 0.120 0.261 0.341 0.176 



Region 1 
Region 2 
Region 3 
Region 4 



Region 1 
Region 2 
Region 3 
Region 4 



